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Abstract — This research developed a breaker depth 
equation based on the characteristic of the potential 
velocity solution of Laplace equation. The breaker length 
equation was obtained using critical wave steepness as 
boundary condition. Wliereas breaker height was obtained 
from breaker height index equation. The equation is in the 
form of linear explicit equation with simple calculation. 
Keywords—breaker depth, critical wave steepness. 

I. INTRODUCTION 

In general, breaker index equation has three shapes, i.e. 
breaker height index —, where H b is breaker height and H 0 

H 0 

is deep water wave height. The second form is breaker 
depth index in the fonn of j^- where h b is breaker depth and 

the third is breaker steepness index in the fonn of ^ where 
L b is breaker length. 

The breaker depth equation developed in this research was 
obtained by performing one of the conservation 
characteristics of potential solution wave of Laplace 
equation, i.e. multiplication between wave number and 
water depth is constant which means that the multiplication 
of wave number at breaker depth is similar with the 
multiplication between wave number and deep water depth. 
Therefore, there is a relation between breaker depth and the 
condition of wave at the deep water. There is a breaker 
length variable at the equation. To eliminate breaker length, 
the critical wave steepness criteria was performed. Using 
this method, an equation is obtained based on analysis and 
the law of conservation. 

II. COMPUTATION OF BREAKER HEIGHT H b 

To perform a computation using breaker index in the form 

of — and —, breaker height H b should be known. The 
Lb hb 

breaker height is obtained using breaker height index 


equation — . There are a lot of breaker height index 

H o 

equations. This research uses a breaker height H b that is the 
average of several breaker height indexes that produces 
adjoining breaker height. The breaker height indexes that 
were used are as follows: 

Komar and Gaughan (1972) 

^=0.56 .,1) 

Singamsetti and Wind (1980) 

— = 0.575m 0 031 (—) - ° 254 .(2) 

H 0 \l 0 J 

Larson and Kraus (1989), 

— = 0.53 (—) °' 24 .(3) 

h 0 \l 0 J 

Smith and Kraus (1990), 

u / u \- 0 . 30 + 0 . 88 m 

— = (0.34 + 2.74m) (—) .(4) 

Hq 'Lo' 

Gourlay (1992), 

— = 0.478 f —) ° 2H .(5) 

h 0 \l 0 J 

Rattana Pitikon and Shibayama (2000) : 

1 

— = (10.02m 3 - 7.46m 2 + 1.32m + 0.55) (—) 5 

H 0 \L 0 J 

.(6) 

At those equations H b is breaker height, H 0 is deep water 
wave height, L 0 is deep water wavelength and m is bottom 
slope.Table (1) presents the result of breaker height 
computation using those 6 (six) equations above and their 
average values. The wave used is the wave with wave 
period T = 8 sec. bottom slope m = 0.005 and water 
depth h 0 = 60 m, whereas deep water wave height H 0 
varies between 0.60 - 1.8 m. 


Table. 1: Breaker height fi-om various breaker height index 
equations 


H 0 

Breaker Height H b (m) 

(m) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Aver. 
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0,6 

0,93 

1,07 

0,91 

1,09 

0,96 

0,93 

0,98 

0,9 

1,29 

1,45 

1,23 

1,48 

1,28 

1,28 

1,34 

1,2 

1,63 

1,8 

1,53 

1,84 

1,57 

1,62 

1,66 

1,5 

1,94 

2,13 

1,81 

2,18 

1,84 

1,93 

1,97 

1,8 

2,25 

2,44 

2,07 

2,5 

2,09 

2,24 

2,26 


Note : (i), number of breaker height index equation 
Table (1) shows that even though each equation provides 
different breaker height, all of them are close enough to the 
average value. 


III. THE COMPUTATION OF BREAKER DEPTH 
h h AND BREAKER LENGTH L b USING THE 
EXISTING EQUATIONS. 

Relation between breaker depth h b and breaker length L b is 
in the shape of implicit equation where there is a 
dependency between those two variables. There is an 
explicit equation to compute breaker depth, i.e. SPM 
equation (1984) and Van Rijn equation (2011). However, at 
Van Rijn there is a parameter that has to be tested, therefore 
this research used SPM equation (1984). The wavelength 
computation cannot be done using dispersion equation of 
linear wave theory, considering at breaker depth the value 
— is quite large, i.e. close to 1, whereas dispersion equation 


of linear wave theory is formulated at very small 



condition. As a consequence, the computation is by doing 
other breaker index equation. 


A. The computation of breaker depth h b and breaker 
length L b with SPM equation (1984) and Miche 
equation (1944) 


Breaker depth index from SPM (1984) is in the fomi of 
explicit equation for breaker depth as an input. 

.(7) 


^6 

H b 


-m 


or h b = 


Hb 




a = 43.75(1 - e“ 19 0m ) 

1.56 

b = 1 + e -19.5m 

H b = breaker height, h b = breaker depth, g = gravity 
acceleration and m = bottom slope. Using (7) h b can be 
calculated using H b as an input, then L b is calculated with 
Miche equation (1944). Whereas breaker index from Miche 
equation belongs to breaker steepness index with the 
following form. 

— =0.142 tanh(—) .(8) 

Lb Lfj 


With the input wave period h b and H b then breaker 
length L b can be calculated using Newton-Rhapson 
iteration method. 


B. The computation of breaker depth h b and breaker 
length L b with Rattana Pitikon et al equation (2003) 
and Miche equation (1944). 

Breaker steepness equation from Rattana Pitikon et al 
(2003) is 

— = (-1.40m 2 + 0.57m + 0.23) (—) ..(9) 

L b \l 0 J 

Using input H b then L b can be calculated explicitly with 

(9), 

l _ Hjl _ 

L b - H 0.35 

(-1.40m 2 + 0.57m + 0.23) (y^) 

Then h b can be calculated using (8), with i Newton- 
Rhapson iteration method. 


IV. THE PROPOSED EQUATION 

The potential velocity of linear wave theory (Dean, 1991) 
is, 0 = Gcoskxcoshkih + z)sinot . This equation was 
obtained by completing Laplace equation with variable 
separation method, where potential velocity is considered as 
consists of Xix) which is only a function of x, Z(z) is only 
a function of z and Tit) is only a function of time, x is 
horizontal axis and z is vertical axis. In (1), Z(z) = 
cosh k(h + z). Applying the velocity potential on small 

. . , i , • 3z(z) a cosh kih+z) 

sloping bottom, and denve -=-= 

r 3x 3x 

\ a/c(/i+z) _ . - . . a k(h+z) _ . 

sinh kyh + z) —--= 0, in this equation —--= 0 in 


ax 


3x 


every z . With k = — and for z = —, the equation can be 

g ( h+ j\ ( h +f) 

. - = 0, which means that -= c, 

dx \ L I L 


written as — 


where c is a constant. From this equation, then if the wave 
moves from water depth h 1 to shallower water depth h 2 


where shoaling took place, then there is a relation, 


h 2 +- 


or. 


h 2 = (hi + yp) — For waves moving 

from deep water with ( h 0 ,L 0 ,H 0 ) to breaker depth with 
ih b , L b ,H b ), there is a relation 



In (10) there are 3 (three) unknowns, i.e. h b ,L b and H b . 
Breaker height can be calculated with breaker height index 
as was done in section 2. Breaker length L b can be 
obtained from the criteria of critical steepness curve. 
Critical steepness from Michel is ^ = 0.142, which states 
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1,2 

1,76 

2,16 

2,17 

2,04 

5,58 

24,36 

38,48 

1,5 

2,08 

2,56 

2,58 

2,44 

6,61 

28,61 

42,17 

1,8 

2,39 

2,94 

2,97 

2,82 

7,59 

32,59 

45,45 


that if - > 0.142 the wave will be breaking, and this 
criteria also applies to breaker point, i.e. breaking occurs 
when — = 0.142. However, in this research, — = - was 

Lb Lb n 

used, this criteria was obtained by studying breaking 
condition with non-linear wave theory, which due to spatial 
limitation could not be shown in this research, and would be 
shown in the next research. Therefore, the breaker length 
becomes, 

L„=nH b (11) 

Substitute to (10), 

i it n b ( i . H 0 ) H b ,. —. 

/l * = —O + yJ-y .(12) 

At the deep water depth, k 0 (h 0 + = a, where 

tanh(cd = 1 , SPM (1984) uses a = n where — = 0.5, in 

L 0 

this research a = 1.17T is used. Substitute L 0 = —to (12) 

k 0 

and a = 1.17T, produces 

h b = (0.55tt - 0.5 )H„ (13) 

With input H b from breaker height index then breaker depth 
h b can be calculated with (13) 

V. COMPARISON OF THE RESULTS OF THE 
THREE METHODS 

The next section will show the result of breaker depth h b 
and breaker length L b computation using the three methods 
mentioned above. The bathymetry data is similar to the 
computation of H b in section 2. The result of the 
computation is as follows. 


Table.2: h b and L b the result of the three methods. 


H 0 

(m) 

H b 

(m) 

K (m) 

M m ) 

(3) 

(1) 

(2) 

(3) 

(1) 

(2) 

Wave Period (X) : 7 sec. 

0,6 

0,92 

1,13 

1,14 

1,07 

2,92 

12,85 

21,59 

0,9 

1,25 

1,54 

1,55 

1,47 

3,98 

17,22 

25,5 

1,2 

1,56 

1,92 

1,94 

1,84 

4,95 

21,15 

28,69 

1,5 

1,85 

2,27 

2,3 

2,2 

5,87 

24,75 

31,45 

1,8 

2,12 

2,61 

2,65 

2,56 

6,74 

28,11 

33,89 

Wave Period (X) : 8 sec. 

0,6 

0,98 

1,21 

1,21 

1,13 

3,12 

13,81 

25,27 

0,9 

1,34 

1,64 

1,65 

1,55 

4,24 

18,55 

29,83 

1,2 

1,66 

2,04 

2,06 

1,95 

5,28 

22,83 

33,57 

1,5 

1,97 

2,42 

2,45 

2,32 

6,25 

26,77 

36,79 

1,8 

2,26 

2,78 

2,82 

2,69 

7,18 

30,45 

39,65 

Wave Period (X) : 9 sec. 

0,6 

1,04 

1,28 

1,28 

1,19 

3,3 

14,68 

28,97 

0,9 

1,41 

1,73 

1,74 

1,63 

4,48 

19,77 

34,2 


Note : 

(1) Couple computation between (7) and (8) 

(2) Couple computation between (8) and (9) 

(3) Computation with (11) and (13) 

The result of h b and L b computation using the three 
methods shows that h b produced by (3) , i.e. with the 
proposed method is quite close with methods (1) and (2). 
Among the three methods, the smallest breaker wave length 
is produced by (3) with quite realistic wave length. 

The computation result of — and — with the three methods 

h b L b 

is presented on table (3). It shows that the value of— from 

hb 

method (3) looks constant against the changes of deep water 

wave height and wave period. At method (1) the value of — 

hb 

is a somewhat affected by the changes in deep water wave 

height but is constant against wave period. The value of — 

hb 

at the third method changes against deep water wave height 

as well as against wave period. The value of — at method 

Lb 

(1) is constant at -, since it was designated. At method (1) 

7T 

the constant is at 0.07, no change against the changes in 
deep water wave height as well as wave period. At method 

(3), the value of — changes against the changes in deep 
Lb 

water wave height as well as wave period. 


Table.3: — and — the result of the three methods 

hb Lb 


H 0 

(m) 

H b 

(m) 

Hb 

h b 

Hb 

Lb 

(3) 

(1) 

(2) 

(3) 

(2) 

(1) 

Wave Period (X) : 8 sec. 

0,6 

0,92 

0,81 

0,81 

0,86 

0,32 

0,07 

0,04 

0,9 

1,25 

0,81 

0,81 

0,86 

0,32 

0,07 

0,05 

1,2 

1,56 

0,81 

0,81 

0,85 

0,32 

0,07 

0,05 

1,5 

1,85 

0,81 

0,8 

0,84 

0,32 

0,07 

0,06 

1,8 

2,12 

0,81 

0,8 

0,83 

0,32 

0,08 

0,06 

Wave Period (T) : 8 sec. 

0,6 

0,98 

0,81 

0,81 

0,87 

0,32 

0,07 

0,04 

0,9 

1,34 

0,81 

0,81 

0,86 

0,32 

0,07 

0,04 

1,2 

1,66 

0,81 

0,81 

0,85 

0,32 

0,07 

0,05 

1,5 

1,97 

0,81 

0,81 

0,85 

0,32 

0,07 

0,05 
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1,8 

2,26 

0,81 

0,8 

0,84 

0,32 

0,07 

0,06 

Wave Period (T) : 9 sec. 

0,6 

1,04 

0,81 

0,81 

0,87 

0,32 

0,07 

0,04 

0,9 

1,41 

0,81 

0,81 

0,87 

0,32 

0,07 

0,04 

1,2 

1,76 

0,81 

0,81 

0,86 

0,32 

0,07 

0,05 

1,5 

2,08 

0,81 

0,81 

0,85 

0,32 

0,07 

0,05 

1,8 

2,39 

0,81 

0,81 

0,85 

0,32 

0,07 

0,05 


Breaker wave length (3) looks very short. This is in 
accordance with the wavelength produced by the 
conservation equation——— = 0. For example, for the 

ignored wave height, relation k 2 h 2 = k^h^ or k 2 = - 1 - 1 

h 2 

applies. For a wave with wave period T = 8 second, moved 
from water depth h ± = 60 m to water depth, h 2 = 2 m, then 

k 2 = = 30 k ± , if k t is calculated using dispersion from 

linear wave theory, it produces k 1 =0.049932 and wave 
length L 1 = 125.835575, produces k 2 = 1.497951 and 
L 2 =4.195 m It shows that equation —— = 0 does 
produce a very short wave length at shallow water. 

From the result above, a simple calculation method of 
breaker parameter can be produced, i.e.: 

a. Breaker height H b is calculated with one of the 
breaker height index equation or by taking average 
values of various breaker height indexes. 

b. Breaker depth h b , is calculated with ^ = 0.81 

c. Breaker length L b j s calculated with ~ = ~ 

’ lb Jr 

VL CONCLUSION 

The result of breaker wave steepness computation with 
breaker wave steepness index equation produces breaker 
wave steepness value that is more or less constant toward 
wave period as well as deep water height. This shows the 
presence of critical steepness wave on a wave curve. 

The proposed equation uses critical wave steepness criteria. 
The equation uses wave condition in deep water and in the 
fomi of explicit equation that is easy to use. In addition, the 
equation is an analytical product based on the law of 
conservation. The critical wave steepness criterion is quite 
important in the development of a simple breaker index; 
therefore, a research is needed on critical steepness wave, in 
the laboratory as well analytical. 
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